Re normalization Group Limit Cycles in Quantum Mechanical Problems 



Erich J. Mueller^ and Tin-Lim Ho^ 

^ Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, New York 14-853 
''Department of Physics, The Ohio State University, Columbus, Ohio 43210 

(Dated: Mar 10, 2004) 

We formulate a renormalizatioii group (RG) for the interaction parameters of the general two- 
body problem and show how a limit cycle emerges in the RG flow if the interaction approaches 
an inverse square law. This limit cycle generates a scaling structure in the energy spectrum. Our 
demonstration is relevant to the Efimov problem in nuclear physics where similar scaling appears. 
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In his famous 1971 paper on the Renormalization 
Group (RG) 1], Kenneth Wilson pointed out that renor- 
malization group equations (which describe how coupling 
constants evolve as one rescales a system) can display all 
of the rich structure of dynamical systems including fixed 
points, chaos, and limit cycles. Yet nearly all examples of 
RG flows in statistical physics contain only fixed points. 
Aside from Huse's 1981 three state chiral model 2j, all 
of the exceptions were found in the last few years. These 
include discrete models of Glazek and Wilson 0| and field 
theory models of LeClair et al. . Amid these revela- 
tions, Albeverio et al. 01 realized that the Efimov effect 
Pi in nuclear physics is associated with a renomalization 
group limit cycle. Here we continue this line of inquiry 
and explore limit cycles in the RG flow of potential scat- 
tering problems (whose origin and structure coincide with 
those of the Efimov problem 6J ) . 

The Efimov effect occurs in a system of three parti- 
cles, when the pairwise interaction is at resonance (i.e. 
the two-body potential is "tuned" to the point where a 
bound state appears just below the continuum). Under 
these circumstances, the three-body system has infinitely 
many bound states, whose energies form a geometric se- 
ries. The ratio between successive binding energies is a 
universal number. Since the deuteron (a bound state of 
a neutron and proton) has such a small binding energy, 
the proton-neutron system is near resonance and the Efi- 
mov problem is relevant to nuclear physics H 0. On 
the other hand, since there are no ways to tune the nu- 
clear potential, exact resonances do not occur in nuclear 
systems, frustrating studies of the Efimov effect. 

The 1996 discovery of Bose-Einstein condensation in 
atomic gases Q has changed this situation. Soon after 
the discovery, Ketterle used a magnetic field to tune 
the atomic scattering to resonance, creating a new system 
for studying the Efimov effect. 

The Efimov problem can be mapped onto a quan- 
tum scattering problem. When the pairwise interac- 
tions are resonant, the important degree of freedom for 
the three-body problem is the hyper- spherical radius R: 
E? = X^i^^j ~ where rj is the position of the 
j'th particle. In this variable, the wavefunction obeys 
a Schrodinger equation with a \/R^ potential Since 



scales as V^, the quantum mechanical problem is 
scale invariant and (as we verify below) can have a RG 
limit cycle. 

In this paper we formulate a renormalization group 
for the interaction constants of the two-particle quantum 
mechanical scattering problem and show how limit cycle 
and fixed point behaviors emerge as the interaction ap- 
proaches and deviates from scale invariance. Although 
two-particle problems are exactly solvable, framing the 
problem in the RG language gives insights which are not 
readily extracted from the exact solutions. For exam- 
ple, our calculation provides a transparent connection 
between scale invariance and limit cycle behavior. We 
believe our simple derivation will be helpful in under- 
standing the generality of the limit cycles in the RG. 

Several recent papers have used RG ideas to discuss 
potential scattering 10, 11] and used these ideas to in- 
vestigate limit cycles |l^ . These authors were interested 
in producing low energy effective field theories (EFT) 
and in constructing short-range regularizations to keep 
energy eigenstates constant. Here, we are interested in 
the global behavior of the RG fiow of the entire set of 
interaction constants as the system is scaled to large and 
short distances. Our method is therefore capable of ex- 
ploring both high and low energy behaviors. 
The quantum mechanical two-body problem: The 
two-body problem in the center of mass frame has the 
form of a single particle scattering problem in a cen- 
tral potential V{r), where r is the distance between the 
particles 0. The energy eigenfunctions ip{Y,E) obey 
the Schrodinger equation {-{h?/M)\I + V{r))tp{Y,E) = 
Eip{r,E), where M/2 is the reduced mass. For our 
discussion, we solely consider s-wave scattering, so -ip 
is only a function of r = |r| and in d dimensions, 
= + {d— l)r~^dr; higher partial waves are treated 
with similar methods. 

The goal of the RG is to analyze a system on differ- 
ent length scale. It is therefore necessary to divide the 
potential into long and short range parts We in- 

troduce a cutoff ro, and let V{r) for r > rg speciiy the 
long distance properties. Short range structure is natu- 
rally encoded in a boundary condition at ro. The most 
general boundary condition consists of specifying the log- 
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arithmic derivative g{r,E) = rip' (r^ E) /tplr, E) at r = tq 
for all E. We express length and energy in dimensionless 
form, X — r/vo (0 < a; < +00), £ — E/e, and U ^V/e, 
where e = It? /Mr^. The same symbol will be used for 
functions of the scaled and original variables [for exam- 
ple, i/'(3^> ^) — ^(^i E), and g{x,£) — xip'{x, £)/'ip{x,£) — 
g{r, E)]. The two-body problem is then specified by 



[~V' + U{xmx,£) 
[xiP'ix,£)/4'{x,£)l^, 



£'ijj{x, £), a; > 1 

5(1, f)- (1) 



Since the choice of Tq is arbitrary, we can choose a dif- 
ferent point, say, Xvo to impose the boundary condition, 
where A is a scale factor. The corresponding equations 
are 



[-V^ + U{x)]ip{x,£) = £'ip{x,£), x>X 
[xij'(x,£)/i;{x,£)]^^^ - g{X,£) 



(2) 



Certainly, the solutions of eq. ^ and Q are identical for 
X > Max(l, A). One can relate g{l,£) to g{X,£) by inte- 
grating the Schrodinger equation, which can be written 
as 

a;^^^ = x^{Uix)^£)-{d-2)g{x,£)^g{x,£f. (3) 

Renormalization: Our rcnormalization group is con- 
structed by recasting the problem in eq. Q into the form 
of eq. 1^ so that the boundary condition is always impose 
at a; = 1. We define scaled variables 



£x = A^f, i)x{x,£x) = '4!{Xx,E), 

g\{x,£x) = g{Xx,£), Ux{x) = X'^U{Xx), 

so that eq.© becomes 

[-V^ + Ux{x)]'ipx{x,£x) = £xipx{x,£x), x>l 
[xip'x{x,£x)/ipx{x,£x)]^^i = gx{l,£\), 



(4) 



(5) 



which has the same form as eq.Q. The "coupling con- 
stants" which parameterize our theory are the potential 
Ux{x) and the boundary condition gx{^,£x), which obey 
"flow equations" 

XdxUxix) = xd,Ux{x) + 2Ux{x) (6) 
Xdxgx{x,£x) = x'{Ux{x)-£x) (7) 
-(rf - 2)gx{x, £x) - gl(x, £x), 

where the derivative of A in dxgxix, £x) acts on both 
subscript A and the argument £x ■ We can classify various 
potentials and boundary conditions by their behavior as 
A increases. 

Flow of Ua(x): To solve eq. 0, we expand Ux{x) in a 
Laurant series 



Ux{x) = Um{X)x"', 

Ti=0,±l,±2,.. 



(8) 



where Um{X = 1) = Um is determined by the initial 
potential U{x). We can explicitly integrate © to find 
Um{X) = t/mA"'+^. We can derive this same result di- 
rectly from the definition of J/a in 

Each variable Um with m 7^ —2, has three fixed points: 
(0, -l-oo, — 00) which are (unstable, stable, stable) for m > 
—2, corresponding to {Um = 0, Um > 0, Um < 0); and are 
(stable, unstable, unstable) for m < —2, corresponding 
to (Um — 0, Um +00, Um ^ — oo). The m — —2 case 
corresponds to the inverse square potential u/x'^ which 
does not flow; i.e. has a fixed line with Um=-2iX) = u. 
Flow of gx- Equation © only involves gx{x,£x) at 
X = 1. Therefore we write gx{£\) = g\i^,£\)- Baring 
pathologies, gx{£x) is a smooth function of £x, and can 
be expanded in a Taylor series 



ffA(fA)=E^=0,l,2,..G.W^A 



(9) 



The function gx{£x) is then characterized by the vec- 
tor G(A) — (Go(A), G'i(A), G2(A), ...), whose components 
obey 

XdxGoiX) - {Y.„^Um{X))-{d~2)Go{X)-Go{\f 
XdxGiiX) = -l-(iGi(A)-2Go(A)Gi(A) (10) 
A9aG,(A) = -(d + 2£-2)Gi{X)-j:,,Gi,iX)Gi-i,{X). 

The last expression holds for £ > 1, and the sum runs 
from ^' = to £. The fixed point solution G* is com- 
pletely specified by the component Gq, 

Gl = -{d + 2Gl)-^ (11) 

" 'd^2 + 2£ + 2Gl ^^^^ 

We can therefore label the behavior of G* by that of Gq, 
which is given by 



Q = {T.mU*m)-{d-2)Gl-G*^ 



*2 



(13) 



Flow of G for non-scale-invariant systems: We 

first consider potentials that do not have a scale in- 
variant component, (i.e. m = —2 term is absent in 
eq.(|Hl). There are then three possibilities for the sum 
J2m^*n = 0, -l-oo, — cx). If the sum is —00, then ifT^ 
has no real solution, and there is no fixed point. We 
will return to this case later. If the sum is -foo, then 
Gq — ±00, (corresponding to a stable/unstable fixed 
point) and G*>q = 0. Physically, this corresponds to 
a strongly repulsive long-range potential. 

For Gq to be finite, we must take the Um = fixed 
point of the potential. The condition that Ux{x) 
as A — > 00 defines a short range potential, and includes 
potentials which vanish beyond a given radius. In this 
case, for d ^ 2 we have 



Gn = and 2 



for d^2, 



(14) 
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FIG. 1: (wide) Cross-sections of renormalization group flows. In each panel, two coupling constants: Um and Go are shown. 
The distinct structures seen when m>— 2,m<— 2 and m = — 2 are illustrated by the three panels. The curves are calculated 
for d = 3 dimensions. Separatrices are shown in bold fl^ . 



with the larger of the two being the stable fixed point. 
(See the left and center panel of figure 1). The fixed 
point vector G* corresponding to each Gq can then be 
obtained from ea. Hll|) to (|12l) '). with stability identical to 
that of their Gq components. As d ^ 2 these two fixed 
points coalesce into one: Gq = 0, which has a saddle- 
point structure. 

To understand the meaning of the fixed point solutions 
Gq , recall the asymptotic form of the wavefunction at low 
energy and large distances |[13 |: 



V'(r) 



\n{r/a) 



2 

d = 2, 



(15) 



where E = fc^/2m and a is the s-wave scattering length. 
Thus, for d 7^ 2, ea. (|15|l implies that for large A, Go = 
{2-d)a'^~^/{{\roY-^+a'^-^). When d > 2 (< 2) and a 7^ 
00, Go flows to (2 — d), corresponding to weak (strong) 
coupling. The unstable fixed point at Go = 2 — d (0) 
is equivalent to a = cx), corresponding to strong (weak) 
coupling. For d = 2, equation (|15|l implies that for large 
A, Go = l/ln(Aro/a). These asymptotics can also be 
determined directly from ea. (|IO|l . 

Flow of G for scale- invariant systems: We now con- 
sider inverse square potentials, U{x) = u/x^. As previ- 
ously shown, C/m=-2(A) has a fixed line solution such that 
Um=-2{X) = u. The first equation of (|1()|) then becomes 

AaAGo(A) = 7.-(d-2)Go(A)-Go(A)2 

= -[Go(A)-A+][Go(A)-A_] (16) 

where A± = {2-d± ^{2 - d)2 + 4u) /2. For u > 

— {d — 2)2/4, the flows resemble those of short ranged 
interactions. There are two fixed points, at Gq = A±, 
with the larger one, A^, being the stable fixed point. 
(See the right panel of figure 1. The interaction C/2(A) 
does not flow as mentioned before). 

For a sufficiently attractive inverse square potential, 
u < —{d — 2)2/4, A± turns into a complex conjugate 
pair, v4_ = 



A+ = ^-r^+i\l\u\-(^-r^\ =a + i(3. (17) 



As a result, ea. (|16|l has no fixed point solution. Instead, 
its solution is a periodic function of ln(A), as seen by 
integrating eq. 116(1 to give 



Go (A) ^a + Pcot (/31nA + 6') 



(18) 



It then follows from ec ljlOl) that G(A) is also periodic in 
ln(A), with the properties 

G(A) = G(AAo), Xo = 2V/3 (19) 

which in turn implies that 



9\{^,£) = .9AA„(l,f). 



(20) 



As A increases, Go(A) goes through ±00. This behavior 
can be rewritten in a more conventional form of a limit 
cycle by performing a conformal transformation which 
maps the real axis onto a unit circle in the complex plane. 



h{\) 



„i<I> 



Go(A)- 
Go(A) - AX 



(21) 



The phase is $ = $0 — 2/31n(A), where, for a pure in- 
verse square potential, tan(2$Q) — /3/(Go(l) — a). The 
new variable h{X) executes simple periodic motion in the 
complex plane. Using the variable the flow of the en- 
tire set of interaction constants is shown schematically in 
figure 2. 

Spectral features: We now show that this limit cy- 
cle leads to a geometric series of bound states and a 
log-periodic structure in the density of scattering states. 
We begin by noting that the eigenenergies {5a ^"■'} of the 
Schrodinger equation in eq. |SJ are related to the original 

eigenvalues {£"(")} in eq.® as = 4"V^^- 

For a sufficiently attractive scale inverse square po- 
tential U{x) = — |u|/x2, equation 1(18(1 holds, so that 
when the system is scaled by A = Ao, the equations of 
motion and boundary condition in (j^J are unchanged. 
Therefore the set of energies {^^l""*} = {A§£*^"-'} must co- 
incide with the original energies {f*^™'}, and each en- 
ergy i?*^"' belongs to a geometric series: for each in- 
teger V, E^'^^/Xl" e In particular, if a sin- 
gle bound state exists, then there must be an infinite 
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FIG. 2: (color) Schematic renormalization group (RG) flows 
near limit cycle. We depict three of the infinitely many dimen- 
sions of coupling constant space. In cyclindrical coordinates 
the {z,p,cj>) axes represent exp((7mi), ^'^v{Um2), and "!>, where 
mi > — 2 and m2 < — 2. These variables, defined in equations 
© and II21II . are respectively relevant, irrelevant, and cyclic. 
The blue horizontal plane and orange vertical cylinder corre- 
spond to the surface Umi = 0, and Um2 = 0. Three RG flows 
are shown: the blue (red) trajectory with Umi = (7^ 0) and 
Um^ 7^ (= 0) spirals in to (away from) the limit cycle; the 
green trajectory lies on the limit cycle with Umi = Um,2 = 0- 

number of bound state whose energies accumulate at 
threshold. Additionally, the density of s-wave states 
Ps{E) = J2i^i^ ~ ^i) must display discrete scale in- 
varience: Ps{\qE) — ps{E)/Xl, so that Ep{E) is peri- 
odic in ln(£^). Coarse-graining over these oscillations, 
Ps{E) ~ !/£' as i? ^ 0. This can be contrasted with 
short-range potentials where Ps{E) ~ I/a/B 16J. 
Strongly attractive long range potentials: We now 
revisit the case where '}2i^ U*^ — > —00 as A 00. Let us 
suppose the Laurant series in ||SJ) terminates at m, and 
let and n = m + 2 > 0. For large A, Go obeys 

AaAGo(A) - -7A" -{d- 2)Go(A) - GoiX)\ (22) 

with 7 = —Ufh{^) > 0. We have neglected terms 
which are small compared to A". As with the case of 
scale invariant potentials, there exists no stationary so- 
lution to these equations and one expects to find pe- 
riodic behavior instead. Although l|22(l is readily in- 
tegrated in terms of Bessel functions, the best insight 
into the behavior of Go comes from noting that Go (A) 
varies more rapidly than A". Thus we perform the con- 
formal transformation in H21() . with j4+(A) — a + i(3 = 
(1 - d/2) + i^-fX"^ - (1 - d/2)2, finding 

Xdxh{X) = -2il3h+[(h^ -l)/2\Xdx\nl3 

= -2i7i/2A"/2/i(A) +C'(A"/2-i). (23) 



The phase therefore behaves as $ = $0 — (4^/7/n)A"/^, 
i.e. it winds around with an increasing rate as A in- 
creases. This changing period is associated with the fact 
that Um{X), continuously increases (see the red curve 
in figure 2), and means that despite the cyclic behavior 
of this is not a limit cycle. 

Ultraviolet properties: If the potential is known to 
arbitrarily short distances, then the scaling equations 
(|l()|l can be integrated to A — 0. Instead of low en- 
ergy physics, one is then exploring high energy physics 
at short distances. The fixed point solutions at short dis- 
tances (which are different from the low energy ones) can 
be obtained similarly and will not be discussed here. We 
only remark that the scale invarient potential U — u/r^ 
has an ultraviolet limit cycle when u < —{d— 2)^/4. 

This work is supported by NASA GRANT-NAG8-1765 
and NSF Grant DMR-0 109255. 

Notes Added: After preparing this paper, we discov- 
ered that the mathematical arguments given here, includ- 
ing the derivation of limit cycle behavior, have previously 
been published by Eugene Kolomeisky and Joseph Stra- 
ley in the context of line-depinning transitions 18] . Due 
to the large overlap between these works, our paper is 
not under consideration for publication in any journal. 
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